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ABSTRACT 

We present near-infrared photometry and spectroscopy, and warm- Spitzer IRAC pho- 
tometry of the young very cool T dwarf Ross 458C, which we have typed as T8.5p. 
By applying the fiducial age constraints lGyr) imposed by the properties of the 
active M dwarf Ross 458A, we have used these data to determine that Ross 458C 
has T c g — 695 ± 60K, log<? = 4.0 — 4.7 and an inferred mass of 5-20 Mj up . We have 
compared fits of the near-infrared spectrum and IRAC photometry to the BT Settl 
and Saumon & Marley model grids, and have found that both sets provide best fits 
that are consistent with our derived properties, whilst the former provide a marginally 
closer match to the data for all scenarios explored here. The main difference between 
the model grids arises in the 4.5 |^m region, where the BT Settl models are able to bet- 
ter predict the flux through the IRAC filter, suggesting that non-equilibrium effects 
on the CO-CO2 ratio are important for shaping the mid-infrared spectra of very cool 
T dwarfs. We have also revisited the issue of dust opacity in the spectra of Ross 458C 
that was raised by Burgasser et al (2010). We have found that the BT Settl models 
which also incorporate a condensate cloud model, provide a better match to the near- 
infrared spectrum of this target than the Saumon & Marley model with f se d — 2, and 
we briefly discuss the influence of condensate clouds on T dwarf spectra. 

Key words: surveys - stars: low-mass, brown dwarfs 



1 INTRODUCTION 

Precise determination of the fundamental properties (i.e. 
mass, age, metallicity and effective temperature, T c a ) of 
isolated field brown dwarfs poses some serious challenges. 
These principally arise because (a) brown dwarfs cool over 
time, and (b) ages are one of the most difficult parameters 
to measure for field stars and brown dwarfs. Together this 
makes breaking the degeneracies between age, mass, lumi- 
nosity and effective temperature for isolated brown dwarfs 
extremely difficult. An object with T e g — 1000K, for exam- 
ple, might be relatively young and low-mass ( e.g. 500 Myr 
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and 20Mj up ), or it could equally be older and higher-mass 
(e.g. 10 Gyr and 60Mj up ). Furthermore, determination of 
Teff for such objects is itself problematic due to the fact that 
the radii of brown d warfs must be estimated using evolu- 
tiona ry models (e.g. iBaraffe et al.|[2003l ; ISaumon fc Marlevl 
120081 ) which to date have yet to be constrained by di- 
rect measurement of field brown dwarf radiQ. The lack of 
age constraints for isolated brown dwarfs makes such es- 
timates particularly uncertain. The influence of metallic- 



1 Although r adii for a handful of brown dwarfs in cataclysmic 
vari able (e.g . Littlcfair ct al. 2007|) and young eclipsing systems 
(e.g. IStassun et alj|2006h have been measured these objects are 
not representative of the field population. 
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ity on the emergent spectra from brown dwarfs further 
complicates matters, and so the calibration of model spec- 
tra using systems with well understood age and metal- 
licity constraints is key to improving our understanding 
of sub-stellar atmospheres. For these reasons, much ef- 
fort has been expended identifying and characterising sys- 
tems for which dynamical masses may be measured (e.g. 
IScholz et al]|2003l ; lDupuv et al.|[2009l ; king et al.ll20ll and 
references the rein) or age and metallicity fiducially con- 
strained (e.g. Pinfield et al.l|2006l; IBurningham et al.ll2009l ; 
iDav- Jones et al.ll2008l.l2010l ; lzhang et alj|201fj )~ 

iGoldman et al l (|2010) report the discovery of one such 
age and metallicity benchmark system, a candidate T8+ 
dwarf ident ified in the UKIRT I nfrared Deep Sky Survey 
(UKIDSS; L awrcncc et al.ll2007l ) that has been shown to 
be a common proper motion companion to the M0.5+M7 
binary system Ross 458AB. Using gy rochronolog y argu - 
ments and the a ge- ac tivity relations of IWest et ah l)2008h , 
IGoldman et"al] |20ld) constrain the age of Ross 458AB 
to less than lGyr. They also highlight its possible mem- 
bership of the Hyades moving gr oup, although confl icting 
radial velocity meas urements by lHawlev et alj J 19971 ) and 
iNidever et all i|2002l) leave this in doubt. However, given the 
spread in age and metallicity seen for members of the Hyades 
moving group, it is not immediately clear that such mem- 
bership would permit any tighter constraints to be placed 
on the pro perties of the Ross 458 system than are possible 
otherwise. iBurgasser et al.l |2010l ) used the Mk /metallicity 
and V — if/metallicity relat i ons of I Johnson fc Appsl (|2009l ) 
and lSchlaufman fc Laughlinl l|2010l ) to estimate metallicities 
of [M/H]= +0.20 ± 0.05 and 0.31 ± 0.05 for the Ross 458 
system. The young (^ 1 Gyr) age of the system, though, 
means that it is not immediately clear that such relations 
can be applied with confidence. We return to this point at 
the end of Section 13.21 

The first estimate of Ross 458C's sp ectral type was 
made through deep methane imaging by IGoldman et al.l 
(2010). It was subsequently and independently identified as 
part of an ongoing large program to identify the coolest 
T dwarfs in the UKIDSS Large Area Survey (LAS) (e.g 
Lodieu et all 120071; iPinfield et all 120081 ; IBurningham et all 



2008, 20091, 12010 j ). As part of that program we have ob- 



tained near-infrared spectroscopy, as well as additional near- 
infrared and warm- Spitzer photometry which we use here to 
estimate the p roperties for this object . 

Recently, IBurgasser et al.l |20ld ) has presented spec- 
troscopy of Ross 458C deriving a spectral type of T8. 
They also compared their spectroscopy to a number of 
oaumon et al l |200rJ , l2007f ) model spectra to derive prop- 
erties. They found that model spectra with significant con- 
densate opacity provided a better match to their data than 
cloudless models. In Section 13.21 we compare our data to 
models from the Saumon & Marley grid, and also to the lat- 
est BT Settl model grid, which also includes a condensate 
cloud model (|Allard et al.ll20ld ). 



2 NEW OBSERVATIONS OF ROSS 458C 
2.1 Near-infrared photometry 

Near-infrared follow-up photometry was obtained usin g the 
Wide Field CAMera (WFCAM; ICasali et al.l 120071 ) on 



UKIRT on the night of 17 th June 2009, and the data 
were processed using the WFCAM science pip eline by the 
Cam bridge Astronomical Surveys Unit (CASU) (|lrwin et al.l 
20041). and archived at the WFCAM Science Archive (WSA; 
Hamblv et al.1120081 ). Observations consisted of a three point 
jitter pattern in the Y and J bands, and five point jitter pat- 
terns in the H and K bands repeated twice. All data were 
acquired with 2x2 microstepping and individual exposures 
of 10 seconds resulting in total integration times of 120 sec- 
onds in Y and J and 400 seconds in H and K. The WFCAM 
filters are on t he Mauna Kea Observ atories (MKO) photo- 
metric system l|Tokunaga et al.ll2002l ). 

The resulting photometry is given in Table [T] It should 
be noted that the follow-up WFCAM data are generally in 
good agreement with the survey data, with the exception of 
the Y band, for which there is a 9a disagreement. The origin 
of this discrepancy is not clear. 



2.2 Warm- Spitzer photometry 

Data were obtained for Ross 458C on 13 th March 2010, via 
the Cycle 6 GO program 60093. Individual frame times were 
30 seconds repeated three times, with a 16-position spiral 
dither pattern, for a total integration time of 24 minutes in 
each of the [3.6] and [4.5] bands. The post-basic-calibrated- 
data (pbcd) mosaics generated by version 18.14.0 of the 
Spitzer pipeline were used to to obtain aperture photometry. 
The photometry was derived using a 7" aperture, and the 
aperture correction was taken from the IRAC handbook- 
The error is estimated by the larger of either the variation 
with sky aperture, or the error implied by the uncertainty 
images. This error is small (Table [T]) and is dominated by 
the 3% error that should be added in quadrature to the 
quoted random errors to account for systematic effects due 
to calibration uncertainties and pipeline dependencies. 



2.3 Near-infrared spectroscopy 

We used JH and H K grisms in the InfraR ed Camera and 
Spectrograph (JRCS; iKobavashi et al.l l2000) on the Subaru 
telescope on Mauna Kea to obtain a R~ 100 JH and HK 
spectra for Ross 458C on 7 th May 2009 and 31 st December 
2009 respectively. The observations were made up of a set of 
four 300s sub-exposures for the JH spectrum and six 240s 
sub-exposures (obtained in an ABBA jitter pattern to facili- 
tate background subtraction), with a slit width of 0.6 arcsec 
delivering a resolution of R~100. In addition we obtained 
a deep H band spectru m using the Near Inf raRed Imager 
and Spectrometer fNIRI; [Hodapp et al.ll2003l ) on the Gem- 
ini North Telescopqj. The NIRI observations were made up 
of twelve 300s sub-exposures observed in an ABBA jitter, 
with a 1 arcsecond slit delivering a resolution of R~500. 

The spectra were extracted using standard IRAF pack- 
ages. The AB pairs were subtracted using generic IRAF 
tools, and median stacked. In the case of IRCS, the data 
were found to be sufficiently uniform in the spatial axis for 
flat-fielding to be neglected. We used a comparison argon arc 



2 http:/ /ssc. spitzer. caltcch.edu/irac/dh/ 

3 under program GN-2009B-Q-62 
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Juki 



Kr 



20.24 ±0.19 17.72 ±0.02 16.69 ± 0.01 17.01 ± 0.04 16.90 ± 0.06 



YwFCAM 



JWFCAM 



HwFCAM KwFCAM 



[3.5] 



[4.5] 



17.54 ±0.02 16.71 ±0.01 17.07 ± 0.03 16.96 ± 0.03 



15.28 ±0.01 13.77 ±0.01 



JWFCAM YwFCAM - JwFCAM JwFCAM ~ HwFCAA 



HwFCAl 



KwFCAM HwFCAM ~ [4.5] [3.5]-[4.5] 



3.53 ±0.19 



0.83 ± 0.02 



-0.36 ±0.03 



0.11 ±0.04 



3.30 ±0.03 



1.15 ± 0.02 



Table 1. Summary of photometry for Ross 458C. The subscripts on the column headings indicated the source of the data, where: 
"UKIDSS" indicates UKIDSS survey data; "SDSS" indicates SDSS DR7 survey photometry; "WFCAM" indicates WFCAM follow-up 
photometry. 



frame to obtain the dispersion solution, which was then ap- 
plied to the pixel coordinates in the dispersion direction on 
the images. The resulting wavelength-calibrated subtracted 
pairs had a low- level of residual sky emission removed by fit- 
ting and subtracting this emission with a set of polynomial 
functions fit to each pixel row perpendicular to the disper- 
sion direction, and considering pixel data on either side of 
the target spectrum only. The spectra were then extracted 
using a linear aperture, and cosmic rays and bad pixels re- 
moved using a sigma-clipping algorithm. 

Telluric correction was achieved by dividing each ex- 
tracted target spectrum by that of an F5 star, observed just 
before or after the target and at a similar airmass. Prior to 
division, hydrogen lines were removed from the standard star 
spectrum by interpolating the stellar continuum. Relative 
flux calibration was then achieved by multiplying through 
by a blackbody spectrum of the appropriate T e g . 

The IRCS JH and HK spectra were joined by scal- 
ing them to match at the H band peak (at 1.58uxn). The 
higher signal-to-noise NIRI H band spectrum was then used 
to replace the H band region of the IRCS spectra, again by 
scaling it to match the IRCS spectrum at 1.58(xm. We note 
here that the NIRI spectrum closely traced the shape of the 
noisier IRCS spectra in this region. 



2.4 Spectral type 

We have assigned a s pectral type for Ross 4 58C following 
th e scheme laid out b y Burga sser et al.l (|2006T l and extended 
by iBurningham et alj ( 20081 ) . Figure [T] shows our combined 
IRCS±NIRI JHK spectrum for Ross 458C along with 



spectra for the T8 and T9 templates 2MASS J04151954- 
093506.6 and ULAS J133553.45+113005.2. Whilst the spec- 
trum of Ross 458C lies roughly between the T8 and T9 
spectra in the J band, it appears to be earlier in type in 
the EbO-if region, and there is considerable excess flux in 
the K band. 

The spectral flux ratios for Ross 458C are given in Ta- 
ble [2] Their values reflect the qualitative comparison with 
the spectral templates described above. The mismatch in 
inferred spectral type between the three indices that are 
not degenerate (Wj, FbO-ff, CYl^-K) makes assigning a 
spectral type for this object problematic. The remaining in- 
dex values are all consistent with a type later than T7, in 
disagreement with the H2O-// and CEU-ftT indices but in 
agreement with the Wj index. We thus allow the Wj index 
to dominate the classification, and assign a type of T8.5p 



(±0.5), reflecting the close agreement between the T8 and 
T9 spectra in the J band, and highlighting the peculiarity 
of the spectrum elsewhere. 

This spectral type is consistent with th at estimated for 
Ross 458C from methane im aging (T8.9, Goldman et al.l 
2010). and also that found bv lBurgasser et al.l J20ld) . How- 
ever, it disagrees with the type of T7 found bv lScholzl (|2010t ) 
using a combination of JHK colour-colour plots and ab- 
solute magnitude arguments. This is driven by the weight 
IScholj l|2010l ') places on the estimates from J — H / J — K 
colour-colour plots, which suggest an early type of T4.0- 
T6.5. Since the colours of Ross 458C are fairly typical of 
other T8± dwarfs discovered to date, it would appear that 
such plots have little utility in spectral typing the latest type 
T dwarfs. 

The disagreement between the Wj index and the 
H2O--H" index that is seen in Ross 458C is broadly con- 
sistent with the unusual spectral m orphology referred 
to as "EbO-ff -early" peculiarity by IBurningham et aU 
|2010ri ). In addition to Ross 458C, two other >T7 dwarfs 
have been identified that display thi s peculiarity: the 
T8p dwarf ULAS J101721.40+01181 7 9 (IBurningham etTal"] 
2008f) and the T7p d warf Gl 229B (|Nakaiima et all Il995l ; 



Burgasser et al.ll2006T ). Several earlier type T dwarfs have 



al so been identified with this peculiarity and are described 
in IBurningham et al.l l^OlObT ). 

It is interesting to note that both Gl 229B and 
ULAS J1017±0118 have inferred ages similar to that for 
Ross 458C, the fo rmer implied by the activity of its M0. 5 
primary star (e.g. I West et al.ll2008T ) and the latter by com- 
parison of near- and mid-inf rared photometry to model pre- 
dictions (|Leggett et alj2010h . However, interpretation of the 
H^O-ff early peculiarity as being directly connected to the 
relative youth of these systems is probably premature, par- 
ticularly in light of the recent observation that amongst the 
ea lier-type objects that displ ay this morphology described 
in IBurningham et all |2010bl ), one has emerged to exhibit 
halo kinematics (Murray et al. in prep). 



3 THE PROPERTIES OF ROSS 458C 

3.1 Luminosity and effective temperature 

We estimate the bolometric flux (Fboi) of Ross 458C fol- 
lowing the method outlined in IBurningham et al.l (|2009T ). by 
combining our JHK spectra (flux calibrated by our WF- 
CAM follow-up photometry) with model spectra (to allow 
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Figure 1. The IRCS+NIRI JHK spectrum for Ross 458C plotted with T8 and T9 spectral standards 2MASS J04151954-093506.6 and 
ULAS JJ133553. 45+113005. 2. Overlaid are the spectral ranges for the numerators (indicated by the index for which they apply) and 
den ominators (D; subscript in dicates the index set) of the spectral typing indices used in Table 1 and outlined in lBurgasser et all (2006) 
and lBurnineham et all d2008l) . 



Index 


Ratio 


Value 


Type 


H 2 0-J 


f 1.165 

/ /(A)eiA 

/•1.28B 

I f(X)dX 

Jl.26 V ' 


0.007 ± 0.004 


>T8 


CH 4 -J 


(■1.34 
Jl.315 /(A)dA 

Jl.26 fWdX 


0.20 ±0.01 


>T8 


Wj 


r-1.23 

f(X)dX 

Jl.18 


0.27 ±0.01 


T8/9 


M.S85 , 
2 f(X)dX 
J 1.26 


K2O-H 


/■1-52 

f(X)dX 

Jl.48 


0.22 ±0.01 


T7 


CH4-H 


/■1.675 

f(X)dX 

Jl.635 
(■160 

Jl.56 /(A)dA 


0.11 ±0.01 


>T8 


NH3-H 


Ji ^3 /(A)dA 


0.71 ±0.01 




CH 4 -K 


f-2.255 

f(X)dX 

J2.215 
<-2.12 

f(X)dX 
J2.O8 k ' 


0.20 ±0.01 


T5 



Table 2. The spectral flux ratios for Ross 458C. The locations of 
the numerators and denominators are indicat ed on Figure [T] The 
NH3 index is not used for assigning a type (see iBurnineham et al.l 
|2008| . and Burningham et al 2010 for a discussion of this), but is 
included for completeness and to permit future comparison with 
other late T dwarfs. 



us to estimate the flux contributions from regions outside 
our near-infrared spectral coverage). We have scaled the 
A < 1.05 urn region of the models to match the flux level 



in our JHK spectrum, whilst we have used the warm- 
Spitzer channel 1 and 2 photometry to scale the 2.5-4.0 
urn and A > 4.0um regions respectively. To provide some 
indication of the systematic effects of our choice of atmo- 
spheric models we have es timated the bolom etric flux us- 
ing both BT Settl models |Allard et al.ll200Sl ) and those of 
ISaumon et ail i|2006l , I2007D . 

We initially selected BT Settl model spectra covering 
the T efi = 500-750K range, logg = 4.0-5.0 and [M/H]= 0.0 
and ±0.3. This range of parameters was selected to ensure 
that we include all likely possibilities for such a young and 
late-type T dwarf, and allowing for the poten tially high 
metallicity suggested by iBurgasser et al. (2010). We then 
took the median as our value for Ftoi , and the scatter in val- 
ues as an estimate of the systematic uncertainty associated 
with our use of theoretical extensions scaled to photome- 
try, amounting to approximately 11%. The uncertainties in 
the photometry used to scale the model spectra contribute 
approximately 4% to the flux error budget. We thus esti- 
mate for Ross 458C that F bot = 6.40 ± 0.80 x 10 _16 Wm~ 2 , 
of which approximately 60% is emitted longward of 2.4um. 
Using the well determined distance to the primary star, and 
incorporating the uncertainty therein, we estimate the total 
luminosity of Ro ss 458C to be 2 59 ± 0.34 x 1O" 6 L . 

We used the lSaumon et all l|2006l . l2007l ) models to fill in 
the missing spectral regions in an identical manner to that 
described above for the BT Settle model grid. We selected 
cloudless models across the same temperature, gravity and 
metallicity range as used for the BT Settl bolometric cor- 
rection, and for simplicity at this stage we have considered 
only a single eddy diffusion coefficient (which parameterises 
non-equilibrium effects due to turbulent mixing - see below) 
of log K zz — 4. This value has been selected as simply a typi- 
cal one for late T dwarfs, and does represent a firm estimate 
for the value of this parameter for Ross 458C, and select- 
ing a different value does not affect our flux estimate, since 
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we scale our model spectra by the observed photometry. In 
Section f3. 21 we find the best fitting value of \ogK zz . 

We find that this choice of models for our bolomet- 
ric correction suggests a value of Fboi = 6.06 ± 0.42 x 
10~ 16 Wm~ 2 , with a 7% contribution to the uncertainty 
arising from the scatter in the model extensions. We find 
that essentially the same contribution to the correction 
arises from the region longward of 2.4um. Our estimate for 
the total luminosity of Ross 458C using the Saumon & Mar- 
ley model set is thus 2.45 ± 0.20 x 10 -6 Lq, which is consis- 
tent with the estimate derived using the BT Settl model 
grid. Both estimates for Lbol are substantially brighter 
than has b een found for the similarly typed late-T dw arf 
Wolf 940B (jBurningham et al.ll2009l ; lLeggett et alj|2010h for 



which 9.79 x 10 7 Lq has been estimated (|Leggett et al.l 



l2010h . 

To determine T e ff for R oss 458C we have used the con- 
straints placed on its age bv lGoldman et al l l|201Ch to esti- 
mate a likely range of radii. As with the bolometric correc- 
tion, we have used two sets of evolutionary models to provide 
some indication of the systematic effects involved in this pro - 
cess. Using predictions from the COND (IBaraffe et alj 2003) 
evolutionary models, for the 0.1 ^ age ^ lGyr range, we find 
that our derived luminosity estimates are consistent with a 
radius in the 0.101-0.120 Rq interval (and log<? = 4.7 and 
4.0 respectively). This results in an inferred T e g range of 
730 ± 25K to 670 ± 20K for our BT Settl based bolomet- 
ric correction, and a range of 720 ± 15 to 660 ± 15 for our 
Saumon & Marley corrected luminosity. 

Using the evolutionary tracks of I Saumon" fc Marlevl 
(2008) suggests a slightly larger radius than the COND 
models, with a radius in the range 0.103-0.124 Rq (and 
log<? = 4.7 and 4.0 respectively), resulting in lower inferred 
values for T eff , with ranges of 720 ± 25K to 660 ± 20K and 
710 ± 15K to 650 ± 15K for our BT Settl and Saumon & 
Marley corrected luminosities respectively. 

Although the BT Settl bolometric correction suggests a 
marginally higher T e g range, the two estimates are consis- 
tent and we adopt T e g = 695 ± 60K as our final estimate, 
with log g — 4.0 — 4.7. For our age constraint of 0.1-lGyr this 
corresponds to a mass range of roughly 5-20 Mj according 
to both sets of evolutionary models. 



3.2 Model comparison 

In Figure [21 we have compared model spectra for a range of 
parameters that straddle our T c a and gravity estimates for 
Ross 458C with the observed spectrum and the mean fluxes 
inferred for the warm- Spitzer channel 1 and 2 photometry. 
On the whole the BT Settl models appear to provide the best 
fits in the near-infrared although, as has been seen before 
(|Burningham et al.| [2009 ) , their K band flux is substantially 
underestimated. Of note also is the difference between the 
two model sets in the 1.57-1.60um region, where missing 
methane opacities are thought to be a significant factor in 
the generally poor fit of model spectra to data in this region. 
Both model sets base their methane opacities on the same 
incomplete line lists, however the BT Settl model grid use 
a statistical estimate of the contributions from hot vibra- 
tional bands in the H and K bands, based on the results of 
iBorvsow etal] (|2003l ), and this may account for the differ- 
ences seen in the H band absorption features. 



The most notable differences between the two model 
sets, however, lie in the 3.5-5um region, where the BT Settl 
models are able to match both the [3.5] and [4.5] fluxes well. 
The Saumon & Marley models, on the other hand, tend to 
somewhat underestimate the flux in the former and overes- 
timate the flux in the latter. Both model sets predict sim- 
ilar absorption strengths for the CO fundamental band at 
4.55um, but a significant difference in the integrated [4.5] 
fluxes arises from a much stronger CO2 band centered at 
4.3um in the BT Settl models than is seen in the Saumon & 
Marley models. Both sets of models use similar CO2 opaci- 
ties and the same CNO abundances, but the BT Settl con- 
sider chemical non-equilibrium effects on the CO2 to CO 
mixing ratio and thus could be expected have higher carbon 
dioxide abundances in the upper atmosphere. 

To examine this possibility in more detail we have recal- 
culated the BT Settl model spectrum for the T c ft = 700K, 
loggr = 4.5, Solar metallicity case under several different as- 
sumptions about CO and CO2 chemistry. The left side of 
Figure [3] compares the 3.5-5.5um region for the BT Settl 
models under the following cases: CE - carbon species in 
full chemical equilibrium; R0 - equal reduction rate for CO2 
and CO; Rl - includes the reaction of CO w ith OH for the 
conve rsion between CO and CO2, as used bv lVisscher et all 
|20ld ) for the Jovian atmosphere; R2 - this is the standard 
BT Settl treatment and is as for Rl but also includes the 
reaction of CO with H2O. The R0 case is roughly equivalent 
to the treatment in the Saumon & Marley models, and the 
right side of Figure [3] shows that this case compares much 
more closely than the two CO2 non-equilibrium cases. This 
supports the assertion that the principal cause of the dif- 
ference between the two model groups in this region is the 
non-equilibrium treatment for CO2 chemistry used by the 
BT Settl models. 

We note that the default BT Settl (R2) case produces 
a between five and six times higher CO2 abundance than 
with the single reaction for converting CO to CO2, and 
30-40 times higher CO2 abundance than in the case where 
both CO and CO2 are quenched at the same level (as in 
the Saumon & Marley models). The increased abundance of 
CO2 does not impact the CO-CH4 chemistry significantly 
since CO remains over 20 times more abundant than CO2 
even in the default BT Settl (R2) model, and CH4 number 
density is similarly unaffected by any of the chemistry mod- 
els since it always remains two orders of magnitude more 
abundant than CO. 

It is clear then that it is appropriate to attribute the 
difference between the models in the 4.3um region to the 
non-equilibrium treatment of CO2 in the BT Settl mod- 
els, and that the strength of the CO2 absorption is also 
strongly dependent on which reaction paths are considered. 
The better fit of the BT Settl models to the observed warm- 
Spitzer photometry suggests that such effects are important 
for determining the emergent spectra of very cool T dwarfs. 
The non-equilibrium treatment for CO2 chemistry is fur- 
ther supported by the recent identification of strong 4.3um 
CO2 absorption in AKARI / IRC spectroscopy of the T 8 
dwarf 2MASS J04151954-0935066 (|Yamamura et all 120101 1. 
although it is not yet established if both reaction pathways 
considered in the BT Settl models are at work in such envi- 
ronments. 

To quantitatively assess the best fitting model we have 
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Figure 2. The JHK spectrum of Ross 458C and mean fluxes inferred from the warm- Spitzer photometry compared to model spectra 
that straddle the properties estimated in Section 13.11 In the case of the Saumon & Marley models, we found that the model grid for 
which logX zz = 6, provided the best match to the data, and this is what has been plotted here. The straight coloured lines indicate 
the mean fluxes of model spectra in the Spitzer photometric bands and plotted to allow comparison with the mean flux from the target 
(straight black lines). 



emplo yed the goodness-of-fit statistic, G, of ICushing et al] 
( 2008) which is defined for a given model, k as 

<=1 

where w is the weight to be assigned to wavelength in- 
terval, i, f is the observed flux with uncertainty a, F is the 
model flux and C is a scaling factor equal to R 2 /D 2 for a 
source of radius R at distance D. We take the value of the 
weight w as the width in wavelength space of each point on 
the spectrum. In all cases we convolved the model spectra 
with the instrumental profile of the observed spectrum of 
Ross 458C such that the resolutions are equivalent. In ap- 
plying this statistic we have treated the total flux captured 
through each warm- Spitzer filter as a single point, weighted 
by the width of the filter. Uncertainties in the value of G 
have been calculated by performing 10000 Monte Carlo sim- 
ulations for each case, randomly offsetting the target spec- 
trum according to the measured uncertainties. 

In Figure [4] we have plotted the values for G for various 
T e ff and logg for both sets of models. We have computed 



C using the known parallax of Ross 458AB and radii con- 
sistent with the values of logg as implied by both sets of 
evolutionary models. We have also calculated G for the case 
where the scaling factor C is unconstrained by parallax and 
evolutionary arguments to determine the T e ff estimate that 
would be found from pure spectral fitting in the case of the 
two sets of models, and the results are plotted in Figure [S] 
We have summarised the best fitting models for each case 
in Table H 

The best fitting model spectra in both the free scal- 
ing fits and the constrained R/D fits suggest properties 
that are consistent with those estimated from our empiri- 
cally determined luminosity for Ross 458C. That the best 
fits are all provided by the Solar metallicity BT Settl mod- 
els appears inconsistent wit h the [M/H] = +0.2 — 1 -0.3 es- 
timates for the system from iBurgasser et al.l (|2010T ). How- 
ever, the Mk /metallicit y an d V — if/metallicity relat i ons o f 
IJohnson fc Appsj (|2009h and lSchlaufman fc Laughhnl (|201Ch 
that were applied in that work are based on calibration from 
a relatively small number of M dwarfs, and Ross 458A lies at 
the early spectral type extreme of the sample. Additionally, 
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Case 


BT Settl best fit 




SM best fit 






T e « 


log 9 


G 


T e « 


log 9 


G 


[M/H] 


= 0.0, fixed C k 


700 K 


4.5 


18.7 + 0.8 


650 K 


4.0 


58.5 + 6.8 


[M/H] = 


-- +0.3, fixed C k 


750 K 


5.0 


30.5 + 2.7 


650 K 


4.0 


38.6 + 4.3 


[M/H] 


= 0.0, free C k 


750 K 


5.0 


16.6 + 0.7 


750 K 


4.0 


33.1 + 2.7 


[M/H] = 


= +0.3, free C k 


700 K 


5.0 


19.8 ± 1.2 


700 K 


4.0 


27.7 + 1.7 



Table 3. Summary of the best fit model spectra and their associated G values. 





Figure 3. Left: A comparison of different treatments of car- 
bon chemistry in the BT Settl model grid for a Solar metallicity, 
T B ff = 700K, logg = 4.5 atmosphere where: CE = carbon species 
in full chemical equilibrium; R0 = equal reduction rate for CO2 
and CO; Rl = includes the reaction of CO with OH for the con- 
version between CO and CO2; R2 = as for Rl but also includes 
the reaction of CO with H2O. Right: A comparison of the R0 
case with the Saumon & Marley model for the same atmospheric 
parameters and logK zz = 6. 



the calibration samples are likely dominated by systems with 
ages typically larger than 3 Gyr. Such systems show little 
evolution in the (V — K s ) — Mk plane with age, and thus 
any spread about the main sequence may be attributed to 
metallicity, however the same can not be said for d warfs 
with the young age of Ross 458. iBurgasser et al.l (2010) also 
noted that the young a ge of Ross 458 may int roduce a bias to 
higher metallicity, and Morales et al.l (|2008l ) found system- 
ati cally higher values o f [Fe/H] for active stars when using 
the lBonfils et all (|2005h relations. 
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Figure 4. The goodncss-of-fit statistic, G, plotted for vari- 
ous model spectra with the value of the scaling factor C fixed 
by evolutionary considerations and the measured parallax to 
Ross 458 AB. For simplicity we only plot the best fitting case 
of log K Z z = 6 for the Saumon & Marley model grid. 



3.3 Dust in the atmosphere of Ross 458C 

In Figure [B] we have plotted the goodness-of-fit statistic, G, 
minimised for the case of a free scale factor, C, for Saumon & 
Marley models in which condensate clouds are present in the 
photosphere and for BT Settl models. The Saumon fc Mar- 
ley atmospheres use the lAckerman fc Marley! l|200ll) cloud 
model, which parameterises the efficiency of sedimentation 
of condensate particles through an efficiency factor f se d- 
Larger values represent faster particle growth and larger 
mean grain sizes. Larger dust grains rain out of the atmo- 



sphere giving rise to optically thin clouds. Smaller values of 
f se( i thus correspond to thicker clouds. 

Since no Saumon & Marley models with dust clouds 
are currently available for log K Z z = 6, we have only fitted 
the near infrared spectrum since a lower value of K zz will 
result in a poor fit at longer wavelengths which would tend 
to dominate the statistic. We found that the best fit for the 
Saumon & Marley models was for a dust sedimentation ef- 
ficiency parameter / se( j = 2, Solar metallicity, T e g = 600K 
and logg = 4.0 (G — 28.03 ± 0.3). These properties are 
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Figure 5. The goodness-of-fit statistic, G, plotted for various 
model spectra with the value of the scaling f actor C evalu- 
ated such that G is minimised in each case (see ICushing et al.l 
l2008t ) .For simplicity we only plot the best fitting case of log K zz = 
6 for the Saumon & Marley model grid. 
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Figure 6. The goodness-of-fit statistic, G, plotted for various 
models, considering just the near-infrared portion of the spectra 
with the value of the scalin g factor C evaluated such that G is 
minimised in each case (see ICushing et~al]|200Sl) . 



Saumon & Marley models (see Figure [2}. We can thus con- 
clude that whilst condensates contribute some near-infrared 
opacity at the low T e g of cool T dwarfs, the use of mod- 
els with clouds thinner in vertical extent than the f se d = 2 
case is preferred. However, it should be noted that the differ- 
ence in the dust opacity between the two models appears to 
be relatively modest. Whilst f se< i = 2 corresponds to very 
dusty atmospheres at the T e g of L dwarfs, it results in a 
roughly similar dust signature to the BT Settl models at 
T cff = 700K. 

We note, also, that the best fitting BT Settl models are 
somewhat warmer than is found when the IRAC photometry 
is included in the fit. This is likely due to the increased im- 
portance of the K band flux in the fit, which is consistently 
underestimated by the cool BT Settl models. 



consistent with what was found by iBurgasser et all l)201fjl l. 
who found that Saumon & Marley models which included 
condensate opacity were able to provide a better fit to their 
near infrared spectra than those without, leading them to 
conclude that condensate clouds are an important opacity 
source in very cool young T dwarf atmospheres. Whilst most 
of the cloud lies beneath the photosphere in this model, up- 
per layers of cloud contribute to the near-infrared opacity, 
with an optical depth of ~ 0.8 in the J band peak, such that 
the cloud is thick enough to define the photosphere in this 
region. 

We have found that the BT Settl model with [M/H] = 
+0.3, T e s = 750K and logg = 4.0 provides a significantly 
better fit (G = 22.05 ± 0.3) than the dusty Saumon & Mar- 
ley models in the restricted near-infrared fits examined in 
this Section. In the BT Settl models the main part of the 
cloud deck also lies well below the photosphere, but the up- 
per cloud layers appear to have a smaller impact than in the 
case of the Saumon & Marley models. This small amount of 
cloud opacity may account for the better match of the BT 
Settl models in the J band peak compared to the cloud-free 



4 COMPARISON TO THE WIDER 
LOW-TEMPERATURE SAMPLE 

To place Ross 458C in context with the wider sample of 
very cool T dwarfs we have reprodu c ed the H — K / H — [4.5] 
colour-colour plot of iLeggett et al.l (|2010T ) in Figure [7] Al- 
though the behaviour of the H — [4.5] colour is likely poorly 
understood (see Section I3.2|l and the models are known to 
produce H — K colours that are too blue, such a plot can 
still provide useful insights to the relative properties of the 
cool dwarf sample. The young age of Ross 458C and its 
location on this colour-colour plot would appear to sup- 
port the assertion that the low-temperature T dwarfs iden- 
tified in UKI DSS to date are d ominated by young, low- 
mass objects l|Leggett et al.ll2~010f ). This is interesting in the 
context of the dearth of late T dwarfs found in the Solar 
neighbourhood compared to what would be expected given 
the mass function seen in young cluster and associations 
|Burningham et al.l [ioiObh . although it is not yet clear if 
this is due to some effect of the birthrate for brown dwarfs, 
an as yet unidentified selection effect or some problem with 
the evolutionary models. 
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Figure 7. Near- to mid-infrared colours of cool T dwarfs com- 
pared to those of the Saumon & Marley model colours, after Fig- 
ure 10 of lLeggett et al.l j2010h . 



5 SUMMARY 

We have presented near-infrared spectroscopy, photometry 
and waxm-Spitzer IRAC photometry of the very cool brown 
dwarf Ross 458C. We assign it the spectral type T8.5p, 
where the 'p' indicates the peculiarity due to an enhanced K 
band flux, presumably reflecting the low-gravity associated 
with such a young object. We have used both BT Settl and 
Saumon & Marley model grids to estimate the flux in regions 
not covered by our near-infrared spectrum and IRAC pho- 
to metry, and by us i ng age constraints placed on Ross 458AB 
bv lGoldman et all (2010) we use the bolometric luminosity 
to estimate that T eS = 695 ± 60K and logp = 4.0 - 4.7. 

We have also fitted various model spectra to the data 
and have found that the Solar metallicity BT Settl model 
grid provide the best match to the data for our derived pa- 
rameter set, although both model sets provide similar es- 
timates for TeH and gravity when used to fit for the pa- 
rameters. The fitted estimates are also consistent with our 
determination from the bolometric flux. 

We note, however, that these similar parameter esti- 
mates are despite significant differences in the predicted 
spectrum in the 4-5 (im region due to differences in ap- 
proach to CO-CO2 chemistry. The BT Settl spectra best 
reproduce the IRAC photometry, suggesting that chemical 
non-equilibrium effects are important in setting the CO-CO2 
ratio, although spectroscopy would be requi red to confirm 
this a ssertion. The FLITECAM instrument (|McLean et all 
2006) on board the Stratospheric Observatory for Infrared 
Astronomy (SOFIA; Becklin 2005) may provide the capabil- 
ity to distinguish the model sets for bright objects such as 
Ross 458C. For fainter objects, this spectral region will be 
restricted to photometric data until the launch of the James 
Webb Space Telescope. 

We have also examined the suggestion by 
iBurgasser et alj |2010l ) that condensate clouds are an 
important source of opacity in the emergent spectrum of 
Ross 458C. We find that the BT Settl models provide a 
better fit to the near infrared spectrum than the dusty 



models favoured by the IBurgasser etafl (|2010l ) fits. Al- 
though in both sets of models the main cloud deck lies well 
below the photosphere, in the Saumon & Marley models the 
upper layers still effectively define the photosphere in near 
infrared spectral regions such as the J band peak, whilst in 
the BT Settl case their influence is much smaller. We thus 

+ .3 

conclude that condensate opacity does indeed appear to 
=4.oaffect the near-infrared spectra of cool T dwarfs, although 
clouds with less vertical extent than is seen for the / se( j = 2 
are preferred. 

=4.5 
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